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Abstract 

A class of linear kinetic Fokker-Planck equations with a non-trivial dif- 
fusion matrix and with periodic boundary conditions in the spatial vari- 
able is considered. After formulating the problem in a geometric setting, 
the question of the rate of convergence to equilibrium is studied within 
the formalism of differential calculus on Riemannian manifolds. Under 
explicit geometric assumptions on the velocity field, the energy function 
and the diffusion matrix, it is shown that global regular solutions converge 
in time to equilibrium with exponential rate. The result is proved by es- 
timating the time derivative of a modified entropy functional, as recently 
proposed by Villani. For spatially homogeneous solutions the assump- 
tions of the main theorem reduce to the curvature bound condition for 
the validity of logarithmic Sobolev inequalities discovered by Bakry and 
Emery. The result applies to the relativistic Fokker-Planck equation in 
the low temperature regime, for which exponential trend to equilibrium 
was previously unknown. 

1 Introduction 



The evolution of many physical or biological systems is characterized by two 
kinds of driving mechanism: diffusion and friction. The competition between 
these two types of dynamics may lead the system to a thermodynamical equilib- 
rium. The purpose of this paper is to study the rate of convergence to equilib- 
rium for a class of linear models that exhibit this kind of behavior. The simplest 
model in this class is the Fokker-Planck equation [23 on the density function of 
the system: 

d tP = Ap + V •(&>), t>0, eet". (1) 
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In this model the density function p depends on the variables (t, £), the dif- 
fusion term is given by Ap and the friction term by V • All sufficiently 
regular solutions of (fTJ) converge in time to a Maxwellian type distribution, with 
exponential rate of convergence [5] . The convergence holds for instance in the 
L^-norm. 

An important generalization of (TTJ, often considered in the mathematical and 
in the physical literature [H H3] , is 

d t p = V -{DiVp + pVE)), (2) 

where D = D(£) is the diffusion matrix and E — E(£) is the energy function 
(E = |£| 2 /2 for {J)). It is assumed that D is positive definite and that 

e- x = / e - E dz 

is bounded. Equation ([2]) admits an unique invariant probability measure, which 
is given by dp, = Poo(Qd£, where p^ — Qe~ E . Non-negative solutions of © 
with unit mass converge to the equilibrium state p^ with exponential rate if 
the matrix D and the function E satisfy an inequality known as the curvature 
bound condition [7] . Let us briefly recall the argument of the proof. In terms 
of ft(i, £) = p{t,£)/pco(£), equation takes the form 

d t h = V ■ (DVh) - DVE ■ V/i , 

or, equivalently, 

d t h = A G h + Qh , (3) 

where A G denotes the Laplace-Beltrami operator associated with the Rieman- 
nian metric G = D^ 1 and Q is the vector field 

Qh = DV log u ■ V/i , u = VdetDe~ E . 

The entropy functional associated with © is given by 

55 [ft] = / hloghdp, (4) 

and satisfies 

|s)[ft]=-g[ft], (5) 

where 

f DVh-Vh , 

m - / — 7 — dp 

is the entropy dissipation functional. Let Ric G and V G denote the Ricci cur- 
vature and the Levi-Civita connection of G. Bakry and Emery proved in [BJ [7] 
that if the curvature bound condition 

Ric G - V G Q > a G , for some a > , (6) 
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is satisfied, then the logarithmic Sobolev inequality T)[f] < (2a) 1 3[f] holds for 
all sufhciently regular probability densities /. Replacing in ([5]) we obtain 



whence the entropy functional decays exponentially as 0(e 2at ). The classical 
Csiszar-Kullback inequality [T3] , 



implies that h converges to 1 as t — > oo in L (d/s) with exponential rate or, 
equivalently, the solution p of ^fy converges to p^ in L 1 (d^) with exponential 
rate. 

Fokker-Planck type equations appear also in kinetic theory and these are the 
subject of the present investigation. Assuming periodic boundary conditions in 
space, the simplest kinetic Fokker-Planck (or Kramers) equation is given by |23j 



d t f+p-V x f = A p f + V p -(pf), t >0, xeT N , per. (8) 



Here / = /(£, x,p) > is the particles distribution in phase-space, with (t, x) 
denoting the space-time variables and p the momentum variable. Equation 
([5]) describes the kinetic motion of a system of particles undergoing stochastic 
collisions with the molecules of a homogeneous fluid in thermal equilibrium. 
From a mathematical point of view, (JSJ is more complicated than ([1]), due to 
the presence of the transport derivative p ■ V x and the fact that the diffusion 
operator A p is degenerate, i.e., it acts only on the momentum variable. Moreover 
equation can be seen as the spatially homogeneous version of © . 

The problem of how fast the solutions of ([5]) converge to equilibrium was solved 
only recently by Herau and Nier [18] and by Villani [26]. Both references es- 
tablish exponential convergence to equilibrium, however by completely different 
methods. In [18] the problem is tackled by spectral analysis techniques for hy- 
poellitpic operators — exponential rate of convergence is implied by the existence 
of a spectral gap in the spectrum of the Fokker-Planck operator — , while the 
proof given in [26] Th. 28] is based on the study of the evolution of a properly 
modified entropy functional. (See [13] for an earlier study of the trend to equi- 
librium for (|8]). In the latter reference the authors prove that convergence to 
equilibrium occurs as fast as Oit^ 1 / 6 ), for any e > 0.) 

In this paper the entropy method is applied to study the trend to equilibrium 
for the following generalization of ([8]): 



where D = D(p) is the diffusion matrix and E — E(p) is the energy function. 
Equation ^ reduces to © in the spatially homogeneous case. The vector field 
v in the transport term is the velocity field; the two most important examples 
for the applications are the classical velocity field 



±J)[h]<-2aQ[h], 




(7) 



dtf + v(p) ■ V x f = V P ( J D(V P / + fV p E)) , 



(9) 



v(p) = p 
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and the relativistic velocity field 



v(p) = — , . 

Setting h — f /e~ E , we may rewrite © in the form 

d t h + v(p)-V x h = A p h + Wh, t > , x G T N , p e R N , (10) 

where A p is the Laplace-Beltrami operator associated with the metric g = D^ 1 

and 

Wh = DVplogu • V p h , u = Vdet D e~ E . 

The main result of this paper is presented in Section [3] It is proved that under 
suitable assumptions on the functions v,D,E , which take the form of geomet- 
ric inequalities involving v,g and W, smooth solutions of (|10p with unit mass 
converge in time to the equilibrium state ft-oo = 1 with exponential rate of con- 
vergence. The proof is based on a generalization of the argument for spatially 
homogeneous solutions outlined before, the main difference being that the en- 
tropy functional Q is replaced by a properly modified entropy functional, as 
proposed recently by Villani [26] . Specifically, the modified entropy used in the 
present paper is given by 

£[h] = kD[h] + a3 pp [h] + 2bJ xp [h] + c3 xx [h] , 

where a, b, c, k are properly chosen positive constants, 

m =11 

is the classical entropy, see (@|, 

3pp[h] = -^p[h] = { [ g(d p h,d p logh)dfidx 
at J t n J r n 

is the entropy dissipation functional (or Fisher information), and 
% x [h] = / / g(A x h,A x logh) dfidx , 

JjN Jt^N 

3xp[h]= / / g(A x h, d p log h)dfi dx , 

where A x h is a modified gradient of h in the x- variable (A x h — d x h when 
v (p) — p)i see Definition OD in Section [3] Note that by Young's inequality, the 
constant a, b, c can be chosen so that £[h] > kT)[h], whence exponential decay 
of the modified entropy entails exponential decay of the classical entropy. The 
reason for choosing the modified entropy in the above form is based on the fol- 
lowing identities, which are proved in Section UTTI Firstly, dD[h]/dt — —3 pp [h]; 



h log h dfi dx 

/ T JV J R N 
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secondly, the mixed derivatives term ^^[/i] appears in the time derivative of 

— 3 P p[7i] = -2Jxp[h] + ... 

where the omitted terms here and below are geometric functions of v, g, W 
and the derivatives of h. Next, differentiating in time the function the 
positive term 3 xa; [7i] shows up: 

^-3 xp [h] = -3 xx [h] + ... 

while the time derivative of 3 XX [h] does not reproduce any of the other func- 
tionals, i.e., 

d ^ 

di Jxx[ 1 " 

The above identities lead to a differential inequality on the modified entropy of 
the form 

p[h]<-S[h} + ... (11) 

provided we assume that the following logarithmic Sobolev inequality in the 
variables (x,p) holds: 

<3*x[h]+3pp[h]. 

In Appendix [B] we show that the latter bound is satisfied when properly defined 
metric and vector field on the product manifold T x R N verify the Bakry- 
Emery curvature bound condition (j6]). The remaining assumptions of the main 
theorem ensure that the omitted terms in the inequality (jllj) give a non-positive 
contribution and Gronwall's inequality completes the proof. The full list of 
assumptions is given in Section [3] They take the form of geometric inequalities 
on the functions g, v, W, which are trivially satisfied for the classical Fokker- 
Planck equation (jHJ. 

It should be noticed that the trend to equilibrium for ([9]) is also studied in [26] . 
The main differences between Villani's approach and the one in this paper are 
the following. In [26] the author exploits the fact that ([9]) can be written in the 
form 

d t h = (A* A + B)h , 

where B = —v(p) ■ V x , A* is the adjoint of A in the Hilbert space L 2 (d^i) and 
A = crV, a = \f~D. The proof of [26j Th. 28] makes crucial use of the iterated 
commutators 

[A,B], [B, [A, B}} , [B,[B,[A,B]]], ... 

in the spirit of Hormander's hypoellipticity theory [TS]. However the use of 
commutators, while natural in the context of regularity theory, presents some 
disadvantages for the problem of convergence to equilibrium. In particular it 
leads to very heavy and sometimes obscure calculations, which, as pointed out by 
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Villani at the beginning of the proof of Lemma 32 in _26 , "might be an indication 
that a more appropriate formalism is still to be found" . The main theorem of 
the present paper is proved using the formalism of differential geometry, which 
helps to clarify the meaning of many long expressions that have to be controlled 
in Villani's work. Moreover, given a diffusion matrix D, a velocity field v and 
an energy function E, one may check directly if our assumptions are verified. 
As opposed to this, one has to find a suitable way to decompose the vector 
fields A,B and the iterated commutators in order to verify the assumptions 
of [ISl Th. 28] . Apart from the different approach to the problem, several ideas 
introduced by Villani in his important work will be adapted (and simplified) to 
the present context, resulting in a less technical proof. 

In Section 3] we apply the main result of the paper to establish exponential 
convergence to equilibrium for the relativistic Fokker-Planck equation [TJ [T51 
116) when the temperature of the surrounding bath is sufficiently small. We 
do not know whether the small temperature assumption is only a technical 
condition or a necessary physical one. However in (2J, see also [5], it is shown 
that spatially homogeneous solutions converge exponential fast to equilibrium 
for all temperatures, although the result is proved not in L 1 but in a weighted 
L 2 norm. 

Finally, we remark that a natural and interesting extension of the results of 
this paper would be to consider the equation on the whole space, x £ M. N , 
with an external confining potential (this is the case studied in [IH1 ) - This 
generalization (and others) will be considered in a separate paper [3], together 
with some more applications of the result. 

2 Set-up 

This section is devoted to introduce the transport-diffusive equation that will 
be the subject of our study, as well as the geometric tools that are needed to 
this purpose. 

Let A^ be a TV-dimensional smooth manifold and M. a smooth M-dimensional 
manifold with a C 2 Riemannian metric g. It is assumed that J\f and A4 are glob- 
ally diffeomorphic to the torus T N and to R M , respectively. This assumption 
will be relaxed in [3] (see also [20] for a formulation of Fokker-Planck equations 
on the cotangent bundle of general manifolds). Let 

h : M x Af x M -> [0, oo) 

satisfy an evolution equation of the following form: 

d t h + Th = A p h + Wh. (12) 

Here h = h(t,x,p), where x — (x 1 ,... ,x N ), p — (p 1 ,... ,p M ) are global coordi- 
nates on Af and M. respectively; dt denotes the partial derivative with respect 
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to t G R, while d x i, d p i denote the partial derivatives in the coordinates (x 1 ,p l ). 
Capital Latin indexes run from 1 to N, small Latin indexes run from 1 to M. 
Denoting X(B), X*(B) the set of smooth vector fields and one form fields on 
a manifold B, then T G X(Af xM), the transport field, whereas W G X{M). 
Finally, A p denotes the Laplace-Beltrami operator on (A4,g). A subscript p is 
attached to differential operators that act on the variables p 1 , . . . , p M only. 

In order to specify the exact form of the fields T, W, some basic facts from 
differential geometry are required. In the following discussion, which is based 
mainly on j22) . we consider only (smooth, time dependent) tensor fields defined 
on M or Af, possibly obtained by projecting tensor fields from Af X M.. We 
also remark that in the rest of the paper we do not distinguish between a tensor 
field defined on N or M and its lift on N x M . Let , P^) denote the frame 
vector fields basis associated with the coordinates x 1 , p % (i.e., X^f — d x if, 

P(i)f = 9 p if, for all smooth functions / on Afx M.) and xi^^P^ their dual one 

form fields; clearly Pu\ and P* are metrically equivalent: g(Pu\,Y) — pJ^(Y), 
for all Y G X(A4). Note that the indexes in round brackets are list indexes 
and not component indexes, that is to say, for each fixed i, P^ is a geometric 
object of the same type (a vector field). Now let ■ ■ ■ v^ N ' denote a set of 
C 3 real valued functions on A4 . We assume that the transport vector field has 
the following form: 

T = v^(p)X (I) . 

We adopt the Einstein summation rule, whereby the sum over repeated indexes 
is understood. 

Remark. When N = M (or more generally when N < M) the functions 
can be thought of as the (non-zero) components of a vector field over M, the 
velocity field. However this interpretation is not necessary and in general not 
very useful, so we will refrain from adopting it. In particular, the use of the 
list index (I) in ?/ 7 ' reminds that this is a scalar function, which affects how 
geometric differential operators act on it. 

For any tensor field R over M, V P R denotes the covariant differential of R, 
where V p is the Levi-Civita connection associated with g (i.e.,V p is symmetric 
and V p g — 0). For a scalar function / on M, V p / is the one form V P /(Y) = 
Y(f), for all Y G X(M). Any vector field Z G X(M) is metrically equivalent to 
the one form field Z* G X*(M) given by Z*(Y) = g(Y, Z), for all Y £ X{M). 
The vector field metrically equivalent to V p / is the gradient of /, which we 
denote d p f: 

g{d p f, Y) - Vp/00 , or V p / - (9 P /)* . 

Using the components gij = g(Pu), P(j)) of the metric in the base P^ (§5 P^ of 
the space of type (2, 0) tensor fields, we may express the action of the Laplace- 
Beltrami operator on scalar functions as 

A p f = ^=d p Jy/\g~\gVd p if) , (13) 
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where g^ is the matrix inverse of gij, i.e., g lk gkj = $ l j and \g\ = detg. There 
is however a more convenient way to express A p /. For this we recall that the 
divergence of a vector field Z G X(M) is the contraction of V p Z, i.e., 

div p Z = V p Z(Pi l) ,P (l) ) . 
We have the well known formula 

A p / = divp (d p f) . 

Moreover by Stokes theorem 

f div p Zy/\g~\dp = 0, (14) 

JR M 

for all vector fields Z G iP(R M , dp). 

Next the definitions of the gradient of a vector field and of the divergence of a 
second order tensor will be recalled. Let V p Z : X*(A4) x X(Ai) — s- R be the 
covariant differential of Z G X(A4). The metrically equivalent type (2, 0) tensor 
field d p Z : X*(M) x £*(.M) -> M given by 

is called the gradient of Z. Moreover, given any type (2, 0) tensor field R, its 
divergence is defined as the contraction of V p i? in the second and third variable, 
i.e., 

div p iZ(F„) = V p i?(n,Pi l) ,P (l) ) 

and thus it is a vector field on M.. The following lemma collects some useful 
identities on the geometric objects defined above. 

Lemma 1. Let /, be smooth real valued functions on A4 and X,Y,Z G 
X{M). Then 

(i) d p {hf 2 ) = hd p f 2 + hd p f 2 ; 
(ii) divp (fZ) = g(d p f, Z) + f div p Z; 

(m) di Vp (d P h ® d P h){Y*) = ag/xtVp/a, n) + A p / 2 a p /i(n) ; 
ff (x,9 p ( 5 (r,z))) = 9 p y(z.,x.) + a p z(n,x t )/ 

(v) g(Y, div p d p 2 /) = 9(Y, d p (A p f j) + Ric(F, p /), 
where Ric denotes the Ricci curvature tensor of g. 

Proof. The proofs of (i)-(iii) are straightforward. The identity (iv) is a conse- 
quence of Koszul's formula applied to the Levi-Civita connection; the proof can 
be found in [22l Ch. 3, Th. 11]. The identity (v) is a direct consequence of the 
definition of the Riemann tensor and is proved for instance in [111 Lemma 1.45]. 

□ 
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We can now define the vector field W. Let E : R M such that 

6- 1 = / e~ E dp 



is bounded. Let d/j, = Oe E dp, a probability measure on R M . Then 

I— e~ E 
d/i = &uy/\g\dp , where u= — -= . (15) 



The triple (A4, g,d/i) is a measure metric space. Let L be the operator in the 
r.h.s. of ([li]). that is 

Lft = A p /i + Wft . 

We require the field W to be such that L is symmetric in the Hilbert space 
L 2 (dn) := L 2 (R M ,d/j,), i.e., 

hLfdfi = / fLhdfi, (16) 

Lemma 2. TTie identity ()16j) is verified if and only ifW = <9 P log it, or equiva- 
lently, W* — V p logu. 

Proof. We have 

/ hLfdfi = Q f h(A p f)uy/\ii\dp+ [ hWfdfl. 

Jm M Jr m Jr m 

In the previous equation we use (i)-(ii) of Lemma [T] to get 

divp (h u d p f) — hu A p f + hg(d p f, d p u) + ug(d p f, d p h) 
and so doing we obtain, by (fT4| . 



hhf dfx = / h(Wf - g{d p f, d p log u)) d/i - / g(5 p /, d/j, . (17) 
Again we use 

divp (fudph) = f u A p h + f g{dph,d p u) + ug{d p f,dph) 
and so we obtain 

/ hLfdp,= / h(W — dp log u)f dfi + / /(A p + dp\ogu)hd/i , 
which implies the claim. □ 
We conclude this section by proving some integration by parts formulas. 
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Lemma 3. The following identities hold true, for all smooth real valued func- 
tions f,h on M. : 



hTfdx = - / fThdx, (18) 

JN JjN 

f hLfd l x = - [ g(d p f,d p h)d!i. (19) 

JR M JR M 

Proof. The proof of (fT8"|) is straightforward. The identity ([!§)) follows by setting 
W = d r hgu in ([TT]). □ 

For the next result we need to recall the definition of inner product of second 
order tensor fields. Given a type (2,0) tensor field R and a type (0,2) tensor 
field S, the inner product R ■ S = S ■ R is defined as 

R ■ S = (R <g> S) (Pi' , P* W , P (i) , P (i) ) . 
Componentwise this means R ■ S — R 1 ^ S^j . 

Lemma 4. For all type (2,0) tensor fields A and Z G X(.M) we have 

g(Z, divp A) dp: = - / A-V p Z*d f i- A(W*,Z*)dfi. (20) 



Proof. Consider the vector field F € 3E(M) defined by F(-) = A(uZ„, •)■ By the 
Leibnitz identity, 

divp Y = g(uZ, div p A) + A ■ V p {uZ*) . 
Replacing in the l.h.s. of (l20l) we obtain 



/ g(Z, divp A)d^ = 6 / g(uZ,d\\ p A)^/\g\dp 

JR M JR M 

= -0 I A-W p (uZ^y^\g~\dp 

Jr m 



/ A ■ (W* ® Z*) dy, - A-V p Z*dn 

jr m Jr m 



which is the claim. □ 



3 Main result 

We begin by stating our assumptions on the functions v^\E and the metric g. 
Let us recall that the Bakry-Emery-Ricci tensor is defined by 



Ric = Ric - VpVK* = Ric - Vl log u , (21) 
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where V?/ denotes the Hessian of / and u is the function (Ti"5)) . As already 
mentioned in the Introduction, Barky and Emery proved in [6) [7] that spatially 
homogeneous solutions of (TT2"j) converge exponentially fast in time to the equi- 
librium state hoc = 1 in the entropic sense (i.e., the entropy functional decays 
exponentially) if the tensor Ric is bounded below by a constant times the met- 
ric g. In the spatially inhomogeneous case, we also need a bound on Ric from 
above. 

Assumption 1. There exist two constants o~2 > ai > such that 

aig{X, X) < Ric(A, X) < a 2 g(X, X) , for all X G X(M) . (22) 

We denote 

(7 = cr 2 - 0-1 > . (23) 

Remark. Anjmportant example is when a = 0. In this case there exists Aei 
such that Ric = Xg. Solutions of the latter equation are called Ricci solitons. 
They play a fundamental role in the analysis of the Ricci flow heat equation, 
see [IT]. 

Before stating the next assumption, it is convenient to give the following defi- 
nition. 

Definition 1. Given a real valued function f on Af x A4, and a point x G Af, 
we denote A x f the vector field over Ai given by 

A x f = (d x i f)d p v^ , evaluated at x G Af. 

The metrically equivalent one form field is given by (A x f)„ = (d x i /)V p v^ . 

We emphasize that A x f G %{M). (More precisely, A x f G X({x} x M) ~ 
X(A4).) Its components in the vector fields basis are given by 

(A x fY =g»d pj v^d x if . 

For the Fokker-Planck equation ([9]) , the manifold M. can be identified with the 
tangent space at all points x G M and A x f coincides with d x f, the gradient in 
x of /. 

Now let us define a symmetric bilinear form A on x R N by 

A(^ v ) = A IJ ^nj, A IJ = g(d p vV\d p v^), Z,r) eR N . 
Note that 

A IJ d x ihd x jh = g(A x h,A x h) . (24) 
Assumption 2. We assume that A is positive definite, 
Mt,0>°, for allO ^ £ eR N . 
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In the next assumption we require the validity of a (weighted) logarithmic 
Sobolev inequality in both variables (x,p). Precisely, define the entropy func- 
tional 



V[h] = hloghdxdfi, (25) 



and 



(26a) 



a r; i f g(A x h,A x h) 
Jxx[h\= I t dxdn, 

JT N xR M 11 

3pp[h]= ( 9(d P h d p h) (2gb) 

Jt n xwl m h 

Assumption 3. We assume that for all smooth functions h : M X M. — > (0, oo), 
such that 

I h dx d[i = 1 , 

JT" xl M 

there exists a > such that the following inequality holds: 

V[h]<^-(j xx [h}+3 pp [h}). (27) 

Logarithmic Sobolev inequalities are extensively studied in the literature and 
several criteria for their validity have been found. In Appendix we give a suf- 
ficient condition for the validity of (|27p . which, in the spirit of our approach, 
takes the form of a geometric inequality on an auxiliary metric defined on the 
product manifold TV x A4. Moreover, as already mentioned in the Introduction, 
the lower bound in (|22|) for the Bakry-Emery-Ricci tensor implies that (|27[) 
holds when h is independent of x (i.e., when 3 XX = 0). 

The previous assumptions suffice if the metric g and the velocity held v are such 
that VpV^ = 0. If this is not the case we need more assumptions, which we 
give after the following definitions. 

Definition 2. Given a real valued function f on M x M., and a point x G TV. 
we denote T> x f the vector field over M given by 

®x/ = (9 x i f)div p dpV^ , evaluated at x € N . 

Definition 3. Given a real valued function f on M x M., and a point x G TV. 
we denote G x f the type (2,0) tensor field over M. given by 

Gxf = {d x i f)dpV (I) , evaluated at x G N '. 

The metrically equivalent type (0,2) tensor field is (C x /)* = (d x i f)VpV^' . 

Next let B, C denote the symmetric bilinear forms on x M. N given by 
B(i,v) = B IJ HiVJ , B IJ = g(dW p d 2 p v^,div p d 2 p v^) , 
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and observe that 

B IJ d x ihd x .,h = g{'B x h,'B x h) , C IJ d x ihd x .,h = t x h ■ (C x /i)* (28) 
Assumption 4. We assume that there exist two constants j3, 7 > such that 
B(£ t Z)<0M<Z,Z), C(t,S)<lMte), forall^eR N . 

The previous assumptions suffice if W lies in the kernel of the Hessian matrix 
of i.e., Vlv^Xw, ■) = 0, for all I = 1, . . . , N. If this is not the case, we 
need a last assumption. Define the vectors . . . K^ N > by 

i.e., componentwise, 

Let R denote the bilinear form on 1SL N x ~R n given by 

R(S, v) = R ij ^ivj , R IJ - g{K (I \ k^) . 

Assumption 5. We assume that there exists a constant w > such that 

R{£,Z)<wA(Z,t), foraineR N . (29) 

Remark. For the Fokker-Planck equation ([§]) we have the following identifica- 
tions: all indexes (small and capital) run from 1 to N — M and 

.'/.,. ».j ■ E = \p\ 2 /2; 

W l = -p l , Bfcy = % , d 2 p v^ = ; 

•A-xf = d x f , T> x f = G x f = , for all functions / ; 

A LJ = Su, B = C = R = i). 

The constants in the assumptions can be chosen as o~\ =02 = a = 1 , /3 = 
7 = lj = 0. A much less obvious example of model satisfying all the previous 
assumptions is given in the next section. 

Before any claim on the asymptotic time behavior of solutions could be made, 
one has to ensure that the Cauchy problem for (fT2|) is globally well-posed. In 
Appendix A we prove a global existence and uniqueness theorem in the L 1 
setting, namely 

0< h G C ([0, 00); L^dxdn)) ; 

we restrict to the case when the dimensions of Af and M. are the same, i.e., 
N = M (which is the most interesting case for the applications) and the metric 
components grow slower than \p\ 2 at infinity. The core of the proof is a general- 
ization of the argument in |17i Ch. 5] , which consists in using the hypoellipticity 
of the Fokker-Planck operator to prove that it is the generator of a dissipative 
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semigroup. We remark that it is possible to prove global well-posedness of the 
Cauchy problem in a much larger class (see the results in [26] for the classical 
Fokker-Planck equation ©), however this is beyond the purpose of this pa- 
per. In the following we assume that the initial datum h[ n belongs to the space 
L 1 C\ L log L(dx d/j) , that 3 X x[hi n } + r 3p P [hi a ] is bounded — see (f2"B"]) — and that h- lTL 
is normalized to a probability distribution: 

\\hn\\Li{dxdii) '■= I h in dxdfi = l. (30) 

JT N xl M 

Our main result is the following. 

Theorem 1. Let the Assumptions HH2I be verified (Assumptions [Zr[7] suffice 
when \7pV^ (W, •) = and Assumptions {J\^ suffice when Vpi/^ = 0) and 
let the initial datum satisfy the aforementioned properties. There exists two 
constants C > 0, A > 0, depending on the parameters o~\, o~2, ui, a, f3, 7, and 
which can be explicitly computed, such that the entropy functional (|25[) satisfies 

V[h](t) < C(3 xx [h in ] + 3 pp [h in ])e- xt . (31) 

Remark. Using ([7]) we have 

\\h - l\\ L i(dxdn) = 0{e- xt/2 ) , as t -> 00 . 

Equivalently, the solution of © converges to the steady state foo ~ e~ E with 
exponential rate in the L 1 norm. 

Remark. Of course there is no loss of generality in restricting to initial data 
that satisfy (|30[) . since (|12[) preserves the L 1 (dxd[i)-iiorm and is invariant by 
the rescaling h —> Mh. Solutions with mass M > converge to the equilibrium 
state hoo = M. 

Remark. It should be emphasized that the entropy method is not suitable to 
obtain the optimal constants C, A for which pip holds and in some cases may 
lead to an exponential rate very far from the real one, see [3S] . The hypoelliptic 
techniques used in [18] are better for this purpose. 

The proof of the Theorem [T] is to be found in Section \5\ In the next section 
we show that the result applies to an important physical model: the relativistic 
Fokker-Planck equation. 



4 Application to the relativistic Fokker-Planck 
equation 

The relativistic kinetic Fokker-Planck equation is obtained from ^ by setting 

p „ I + p <8> p 



° lp) = 7T7W D =vTf#' E = ° VTTW ' (32) 
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We restrict to the three dimensional problem: ieT 3 ,peM 3 (thus capital and 
small Latin indexes run both from 1 to 3 in this section). v(p) is the relativistic 
velocity, D is the relativistic diffusion matrix and y/l + \p\ 2 is the relativistic 
energy. We set the rest mass of the particles and the speed of light equal to one. 
The equilibrium state is given by the Jiittncr distribution 

Jeip) = Ze- e V^^ , (33) 

where Z is a constant (fixed by the mass of the system) and 9 is a positive 
parameter which, up to a dimensional constant, coincides with 1/T, where T 
is the temperature of the surrounding bath in which the particles are moving. 
Although the interest in the relativistic Fokker-Planck equation has increased 
substantialy in recent years [TJ [TSl [H] , the relativistic theory of Brownian mo- 
tions is an old classical topic [14]. In this section we prove that solutions of 
the relativistic Fokker-Planck equation converge with exponential rate to the 
Jiittncr equilibrium, provided the parameter 9 is sufficiently large, i.e., for a 
sufficiently small temperature. We do so by showing that all the assumptions 
of Theorem Q] are satisfied. To this purpose we first normalize the solution by 
introducing h — fe E and rewrite the relativistic Fokker-Planck equation in the 
form (|12p , where the metric g is given by 

9ij = Po($ij - ^ir) , po = V 1 + IpI 2 ■ 

Po 

Note that g is conformal to the hyperbolic metric 5y — pipj/p^ and that the 
matrix inverse of g^j is 

9 ij = ~(S ij +pV), 
Po 

where the indexes of the variables pi are raised and lowered with the Euclidean 
matrix, e.g., pi — b^pK Since det g — p , the function u = e~ E / y/det g is given 

by 

e -epo 



u = 




We begin by showing the validity of Assumptions[2]and[H which are independent 
of the energy function E. A straightforward calculation shows that the bilinear 
form A IJ is given by 

A IJ = g(d p vW,d p vW) = fidpivWd^vW = \(S IJ - P -^-) . 

Po Po 

Thus Assumption [2] holds, because 

mu) = - ^#) > iei 2 (i - K) - K > o , 04) 

Po Po Po Po Po 
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for all ^ i £ M 3 . The bilinear form B IJ and C IJ are given bjQ 

u = 496pg - 9030pg + 1035pg - 2b „ 
16p 3 

25 - 1035pg + 10551^ + 1610pg + 729p§ J7 

Since, by (jMj) . 

l£| 2 <i^(£>0 (35) 



holds, we have 



where PsCpo) is a polynomial of degree 8. Thus it is clear that there exists f3 > 
such that B(£, £) < £) and by the same argument, there exists 7 > such 

that C(£,£) < jA(£, £). We conclude that Assumption |4] is satisfied as well. As 
to Assumption [SJ the bilinear form R IJ is given by 

i?/J = (1 + 2%^ _ „ + _ + 

and arguing as before it is easy to prove that £) < £), for some 

positive w. Let us now take care of Assumption [TJ The Ricci tensor of g and 
the Hessian of log u are given by 



p . 1 4 + 15pg 

Ric y = T3 (3dij Wj) 

4 Po 



(Vj logu)y = — s ((4 + 49 Po )6 l3 ^ ^ ^- 9ij ) . 

^Po Po 

Therefore the Bakry-Emery-Ricci tensor reads 

gr 1 , M , , fl u 6^-12^ + 2^-1 

Rlc *i = 772 (-(! + 4 %)% + 9ij) ■ 

^Po Po 

It is straightforward that the bound from above in (|22|) is satisfied, for all 9 > 0. 
However the lower bound is satisfied if and only if 9 is sufficiently large. To see 
this we use that 

\X\ 2 < p Q g{X,X) , for all X € K 3 , (36) 



whence 



4Po 



The remaining calculations in this section have been carried out with MATHEMATICA. 
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and the minimum of the function on square brackets is strictly positive if and 
only if 9 is sufficiently large (e.g., 9 > 4 suffices). It remains to check the validity 
of the logarithmic Sobolev inequality (1271) , where in this case 

3pp[h] = I — (\V p h\ 2 + \p-S7 p h\ 2 )dxdfi . 

To prove (j2"?| we use Theorem 0] in Appendix B. We introduce the metric 

G = Gudx 1 dx J + g i jdp l dp j , 
where Gij = Au — p\{^ij + pipj) is the inverse matrix of A IJ , the function 

u= u 



and check that there exists a constant a > such that 



p 

The Ricci tensor of G is given by 



Ric G (Z, Z) - Vp log U(Z, Z) > aG(Z, Z) . (37) 



Ricfj = 


13 a . 


19^-7^ 




pI 


Ricfj = 


o, 




Ric? = 


3 r 

2p 2 o ij 


2hp 2 - 3 
2pl 9 «- 



The Hessian of log U is 

2 23 + 29po 2 5pg-2 
(V log U)u = (2p Q 5u 3 — 

4 Pq 

(V 2 logU) iJ = 0, 



fu2| rn 23 + 9 Pa 6flpg + m P j + 29 Po + 69 

(V logt0tf = — 2 % 7-3 fftf- 



Using the previous formulas, the bound (|36p and 

|X| 2 < GljX ^ XJ , fo r all lei 3 , (38) 

it is straightforward to prove that the curvature bound (|37p , and thus the loga- 
rithmic Sobolev inequality (|27p . holds when the constant is sufficiently large. 
In conclusion all the assumptions of Theorem Q] are satisfied for the relativistic 
Fokker-Planck equation provided the constant 9 is large enough, or, equiva- 
lently, the temperature of the surrounding bath is sufficiently small. Therefore 
Theorem Q] yields the following result. 
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Theorem 2. Let < f m be an initial datum of mass M > for the the 
relativistic Fokker-Planck equation, i.e., for ([9]) with Q32[) substituted in. Denote 
by 3b,m the Jiittner distribution (|33p with mass M . Then there exists 8q > 
such that for all 8 > 8q there exists two positive constants C, X, depending on 8, 
such that the solution f of the relativistic Fokker-Planck equation satisfies 



Remark. In [5] the precise values of the constants C, A, 8q are found in the case 
of spatially homogeneous (SH) solutions. Moreover it is shown that, at least 
within the class of SH solutions, exponential convergence to equilibrium holds 
for all 8 > in a suitable weighted L 2 norm. 

5 Proof of the main result 

In the rest of the paper the following abbreviations will be used: 



Moreover the measure dx d[i will be omitted in the proofs. 

In the following we assume that ft is a positive smooth solution of p^|) . The 
proof of the result for L 1 non-negative solutions with finite entropy is obtained 
by a standard approximation argument, see [H [10] for examples of this proce- 
dure. 

Recall that 



\\f - JeMv < Ce 




■ ■ ■ dx d/i 



and 



h = log h . 




and define the mixed derivatives term 




Given four constants a, b, c, k > 0, we define the modified entropy as 



£[h] = kD[h] +aJ pp [h] + 2bJ xp [h] + c3 xx [h] . 



We divide the proof in three subsections. 
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5.1 Evolution of the modified Entropy 

Our first goal is to study the time evolution of the modified entropy, by com- 
puting the time derivative of D, 3 XX , 3 xp and 3 PP . 



Lemma 5. The following holds: 



j t ®[h] = -3 pp [h] 



Proof. We compute 

-£8[h] = J d t h(l + h) = - J (l + h)Th + J Lhdxd f i + J hLh. 
By (fT51) and (JTHJ), the hrst two terms vanish and 

J hLh = — J g(d p h, d p h) . 



Lemma 6. The following holds 

d_ 

dt 



3 pp [h] = — 2J a: j ) [/i] — 2 J h Ric(d p h, d p h) dx d[i 
- 2 J hd%h- \7 2 p hdxd[i . 
where Ric is the Bakry- Emery- Ricci tensor (|21[) . 



Proof. We compute 
d 



dt ^pp[ h ] = 2 / g(d p h,d p d t h) - J g(d p h,d p h)d t h 

= -2 { g(d p h,d p {Th))+ [ g(d p h,d p h)Th 



<? 
+2 / g(d p h,dp(Lh)) ~ / g(d p h,d p h)Lh 



□ 



We claim that V + (} = —23 xp [h]. We prove this using the coordinates repre- 
sentation. From one hand 

g(d p h, d p {Th)) = g l: >d p ,h(d p] v {I) )d x ih + g ij d r h v {1) d p] d x i h ; 

on the other hand, integrating by parts in the x variable, 

= -2 / g l3 d p] hd p id x ihv {1} = 2 f g ij d p] d x i h v {1) d p ih . 
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Thus 

9 + = -2 / g l3 d p ,h{d p] v^)d x ih = -2 / g(A x h,d p h) . (39) 



The term X is 

* = 2 y 0(9^, 9 P (A P ^)) + 2 y g(d p h, d p {Wh)) = *, + * 2 . 
By (v) of Lemma [T] and (|2ll|) we have 

*i = - 2 y Ric^/i, 9 p ft) + 2 y g(d p h, dw p d%h) 

= - 2 y Ric(d p 7i, -2 J d 2 p h-V 2 p h~2 J d 2 h(W*,V p h) . 
Moreover by (iv) of Lemma [TJ 

* 2 = 2 J g(d p h,d p (g(d p h,W))) = 2 J d 2 h(W^V p h) + 2 J d p W(V p h,V p h) 
Summing up and using the identity 



d 2 h = h d 2 h + d p h ® d p h 



we obtain 



4 = -2 J hRic(d p h,d p h) - 2 J hd 2 h- V 2 p h -2 J V 2 p h(d p h,d p h) . (40) 
Finally, by and (iv) of Lemma [TJ 

* = [g(d p h,d p (g(d p h,d p h)))=2j%h(V p h,V p h), (41) 
which cancels out with the last term of (1401) . The claim follows summing 

up gehuj. □ 

Lemma 7. The following holds: 

-jpx P [h] = -3xx[h]- / hRic(A x h,d p h)dxdfi 

~?,J hd 2 h ■ V p (A x h)* dx dfi + J g(d p h, r B x h) dx dfi 
+ 2 J hd p h- (e x h)*dxdn + J e x /i(W*,V p h) dxd^i. 
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Proof. We have 

^px P [h] = f g(A x d t h,d p h) - J g(A x h,d p h)d t h + J g(A x h,d p d t h)) 

r g(A x h,d p {Th))- I g(A x (Th),d p h)+ [ g(A x h,d p h)Th 



<?1 ^2 <?3 

g{A x h,d p h)Lh+ [ g(A x h,d p (Lh))+ [ g(A x (Lh),d p h) 



* * 

Now we claim that 

V = Vi+V2 + Vs = -3**[h]. ( 42 ) 

In fact, using the coordinates representation the first term of 's? can be rewritten 
as 

9i = - J g ij d p iV (I) d x ihd pi v w d x jh- J g ll d p ,v {I) d x ihv {J) d p] d x .,h 

It is clear that ^?±a — —3 xx - Integrating by parts in the x variable we obtain 

<?m = J g rj d p iv {1) d x id x .,hv {j) d p] h. 

In the previous expression we use the identity 

d x id x .ih = h~ 1 d x id x jh — h~ 2 d x ihd x jh 
and so doing we obtain 

g(A x (Th),d p h) - [ g(A x h,d p h)Th = -<? 2 - <? 3 . 



This proves (|42|) . It remains to study the integrals 0, We begin by 

applying (fT9|) and (iv) of Lemma [1] to <0>: 



= / g(d p h,d p (g(A x h,d p h))) = J d p (A x h)(V p h,V p h)+J d 2 p h{{A x h)*,V p h) . 

(43) 

As to we first split it as 

* = J g(A x h,d p (A p h)) + J g{A x h,d p (Wh)) = *! + * 2 . 
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By (v) of Lemma [T] and (|2T)|) we have 

g{A x h, div p dph) - J Bic(A x h, d p h) 
dlh-V p {A x h)* - [ d 2 h{W,,{A x h),)- I Bic(A x h,d p h) 



Likewise 



* 2 = / g(A x h,d p (g(d p h,W))) = J d 2 h(W*,{A x h) ir ) + J V p W*(d p h,A x h) . 
Summing up, 

* = - / ^c(A x h, d p h) - J dpi- V p {A x h), . 

In the second integral we replace 

d p h = h dph + dph (g) d p h 

and we get 

* = -J hWc(A x h,d p h)- J hd 2 h-V p {A x h),~ J dp{A x h){Vph,Vph). (44) 

Note the the last term in the r.h.s. of (|4"4")) cancels out with the first term in 
the r.h.s. of (|4U)) . We now work out the term 4k. First we rewrite it as 

* = J L{d x ih)g{dpv^\dph) . 

Then by (TIT?)) and (iv) of Lemma [1] we have 

* = - f fig« W (V P h, V P d x i h)- J dfh(V p vW , V p d xI h) = *i + *2. (45a) 

In 4b i we apply the identity 

Vp/i <8> V p d x ih = Vp{V p hd x ih) — \7 2 hd x ih 
and (l20l) to obtain 



A = - / d 2 v {I \\7 p h,\7pd x ih) = - J d 2 v {1) -\J p h®\J p d x ih 
d 2 pV {I) d x ,h-V 2 ph- [ d 2 p v {1) -VpiVphd^h) 



= / ft figft • (C x ft)* + / g{dph,V x h) + / % v W(W„d t ihV p h) . (45b) 
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Note that the last term in (|45b|) can be rewritten as 

J d*vW(W*,d x ihV p h)= f e x h{W*,V p h) . (45c) 

In 4*2 we apply the identity 

V p d x ih = hV p d x ih + d x ih\J p h 

to obtain 

* 2 = - J d%h(V p vW,d x ihV p h)- J hd 2 p h{V p v^\v p d x ,h) 

= - J d 2 h((A x h)*,S7 p h)- J hd 2 h- \J p v {1) ®\7 p d x ih = * 2A + *2B ■ 

(45d) 



Note that 4f*2A cancels out with the second term in the r.h.s. of (l43l) . In 4tk 2 B 

we use 

V p v {1) ® V p d x ih = V p (A x h)* - Vlv^d x ih 

to finally obtain 

+2B = - J h d 2 p h ■ V p (A x h) t + J hd 2 h- (G x h) t . (45e) 

The claim follows by (|32" ]) -(|4l) j) . □ 
Lemma 8. The following holds: 

= — 2 J hd p {A x h) ■V p (A x h)*dxd[j,-\-2 J g(A x h,15 x h) dx d/i 

+ 4/ Ul (4S).(a).^ + 2/w»'.,w).)**. 

Proof. The proof is very similar to that of Lemma [7] First we compute 
d 



df -3xx[/i] =2 / g(A x d t h,A x h) - J g{A x h,A x h)d t h 

= -2 I g(A x (Th),A x h)+ { Thg{A x h,A x h) 



v o 
Lhg(A x h,A x h)+2 [ g{A x (Lh),A x h) 



We claim that 



<? + = . (46) 



2:-S 



In fact, by (15)1 

= -J hT(g{A x h,A x h)) = -2 J hg(A x (Th) , AJi) 

= -2 J g(A x (Th),A x h) + 2 J ' (Th) g(A x h,A x h) => = -<? . 
By (fl~9f and (iv) of Lemma Q] the term «k can be rewritten as 

* = 2 [ d p (A x h)((A x h)*,V p h) . (47) 

Likewise, 

4 = 2^ L(d x rh)g(d p v^,A x h) = - 2 / g(d P h,d p (g(d p v^ ,A x (d x ih)))) 

= -2 J dy i \{A x h)^V p d x ih)-2 J d p (A x h){V p v^\V p d x ih) 

= *i + * 2 ■ (48a) 

Using the identity 

(A x h)* ® \7 p d x ih — S7 p ((A x h)*d x ih) — d x ihV p (A x h)* , 

we may rewrite 4k i as 

*i = - 2 J d p v (I) {{A x h)«,V p d xl h) = -2 J d 2 v^ ■ (A x h)*®V p d x ih 
= -2 J d 2 v^ ■ V p ((A x h)*d xI h) + 2 J d 2 p v^ ■ V p (A x h)*d x ih . 
Applying (|2U| to the first term in the last line we get 

4>i=2 J g(A x h,T> x h) + 2 J d 2 vW(W*,(A x h)*d x ih) + 2 J hd p (A x h) ■ (Q x h)* . 

(48b) 

The second integral in the right hand side can be rewritten as 



2 J d 2 p v {I \W^ (A x h)*d x ih) = 2 J e x (W*, (Ah)*) . (48c) 
Using the identity 

V p d x ih = hV p d x ih + d x ih\J p h 

the term 4ko becomes 

* 2 = -2 J hd p (A x h)(V p v ( - I \V p d x ih)-2 J d p (A x h)(V p v^,d x ihV p h) 
= **2A + *»2B ■ (48d) 
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Note that 4*2 b cancels out with X In 4*2 A we use 

d p (A x h)(VpV W ,V p d x ih) = d p (A x h) ■ V p v {1) ® V p (9 x /ft 

and 

V p w (/) ® V p d x ih = V p {A x h)» - V 2 p v {I) d x ih 

to obtain 

=-2 /" hd p (A x h)- V p (A x h*) + 2 J hd p (A x h) ■ (C x /i)* . (48e) 
Summing up (|4l))) - (|4"5)) concludes the proof. □ 

5.2 A differential inequality for the modified entropy 

Recall that 

£[h] = fcJD[ft] + al pp [h] + 2b3 xp [h] + cJ xx [h] . 

In this section we prove that, under suitable conditions on the constants a, 6, c, k, 
the modified entropy satisfies 



£[h] >kD[h], —£[h] < -d(a3 pp [h] + 2b% p [h] + c3 xx [h}) , (49) 



where d is a positive constant. In particular, the first bound shows that expo- 
nential decay of the modified entropy implies exponential decay of the entropy. 

The bound from below is easily established. 
Lemma 9. Assume b < \[ac. Then £[h] > kD[h]. 

Proof. By Young's inequality, for all e > we have 

g(A x h 7 d p h) > ~eg(d p h,d p h) - ^-g(A x h 7 A x h) , 

whence 2b3 xp > —2be2 pp — (b/2e)3 xx and so 

£[h] > kD[h] + (c - 3 xx [h] + (a - 2eb)3 pp [h] > kD[h] , 
provided b/c < 2e < a/b. □ 

The bound from above, which requires the assumptions of the main theorem 
(except Assumption [3]) , is more complicated. Since 

a3 pp [h] + 2b3 xp [h] + c3 pp [h] < max(a + b, b + c)(3 pp [h] + 3 xx [h]) , 

it suffices to prove the following. 
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Proposition 1. Let the Assumptions [7J [H [7) [3 hold. There exists a region 
Q, C R 3 such that fl C {(a, b,c) : b < ^/ac} and, for all (a, 6, c) <E ri, £/iere exists 
d > smc/i £/ia£ 

< -d(3f ra [/i] + 3 pp [/i]) . (50) 

Remark. The best constant in the inequality (|50|) may be written as d = sup^ d, 
where f2 is the largest region for which Proposition [T] holds. We refrain from 
computing it explicitly, since the method we use is anyway unsuitable to obtain 
the optimal rate of decay of the entropy. 

The proof of the proposition is based on the following lemma. 
Lemma 10. For all constants e±, . . . , e±o > we have 

j^ppVA < 2e 1 3 xx [h] + - 2a^j 3 pp [h] - 2Q 2 p 



(51) 



e 2 cr + e 3 (Ti + 2e 5 7 + e 7 ui + e 6 (3 - 1 3 xx [h] 



£2 £3 



£6 



£7 



2e 4 



2 

PP 



2 

xp 



(52) 



^xx[h] < ( 4£ 87 + ^- + 2£ 9 /3 + 2e 10 oj + — !- ) 3 xx [h] + ( — - 2 ) Q 2 

at \ Z£g ^£10/ \£8 / 



(53) 



where 



Q 2 p = / hd 2 h-V 2 hdxdfi , Q 2 xp = / hd p {A x h) ■ V P {A X K) dxdfi 



Proof. The inequality (|51j) is a straightforward consequence of Lemma [6j the 
inequality 3z P [/i] > —Si3 xx {h} — (4ei) _1 3 pp [/i], and Assumption [TJ We now 
prove (|52"|) . Using the identity 



Ric(A. x h, d p h) =Ric ( \[e2A x h H — d p h, ^/e^,A x h H — d p h, 



/4e 2 



/4e 2 



— £2Ric(7l :! ;/i,yi s /i) — - — Kic(d p h,d p h) , 

4^2 

together with Assumption [T] and 3a. p [/i] > — e^3 xx [h] — (4e 3 )~ 1 3 pp [/i], we get 



- J h Ric(A x h, d p h) < (e 2 o- + £ 3 (Ti 
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By Young's inequality 

-2 f hd%h-V p (A x h)* < 2e A J hdfi-Vlh+^-J hd p (A x h)-V p (A x h)*. (55) 
By Young's inequality, ([28]). Assumption [4] and (|24l) 

2 J hd 2 p h- {e x h)* <2e 5 J h e x h ■ (e x h)* + ^.f hd fi- v ^ 

<2e 5l 3 xx [h] + ^-Jhd 2 p h-V 2 p h. (56) 

Likewise 

/ gidph^h) < e 6 / hg^Ji^Ji) + — 3 pp [h] < e 6 pj xx [h] + -?-3 PP [h] . 
J J 4e 6 4e 6 

(57) 

Finally by Assumption [SJ 

<Z x h{W*,V p h) = y g(K^d x ih,d p h) 

<e 7 J hg{K^\K^)d x ,hd x .,h + J g(d p h,d p h) 
< e 7 w3 xx [h] + - — 3pp[h] . (58) 

4£7 

Using the inequalities (f5~4" ]) - (|5"8")) in Lemma [7] concludes the proof of (|5"2")) . The 
proof of (|53|) is similar. Reasoning as before one can prove that 

4 / hd p {A x h) ■ (Q x h)* < 4e sl 3 xx [h] + — ; 

2 /" g(A x h,-B x h) < (2e 9 (3 + ^)3 xx [h] , 

2 y Q x h{W*, (A x h)*) = 2 y hg(K^d x ih,(A x h)) < (2e w u; + ^-)3 xx [h] 

and substituting in Lemma |8] completes the proof. □ 

Remark. We are going to apply Lemma [10] for special values of the constants 
£i, . . . , eio- In its generality, Lemma HU1 could be useful to improve the constant 
d in (E 



Proof of Proposition^ In the inequalities (|51[) - (153l) we set 

£1 = (2a) _1 , £ 2 = £ 3 = £6 = £7 = t(o"2 + /3 + , 
g 

£ 4 = -(2 + /3 + 167 + w) , £ 5 = (87)-\ £ 8 = 4,£ 9 =£ 10 = 
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So doing we obtain 



—3pp[h] < a 1 3 xx [h] + (a - 2<7 1 )l pp [h] - 2Q 2 pp , 



PPvi -~**pp 

16s 
~7~ 
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i>ii 



—0 2 



dt 
where 

Therefore 
d 



-ir^xxih] < s3 xx [h] - -( 



si = o- 2 + + uj , s 2 = 2 + cr 2 , s = 2 + f3 + I67 + u 



a^3 pp [h] + 2b^3 xp [h] + c-^3 xx [h] 



< [—k + a(a — 2<Ti) + 2bs\S2\ 3 pp [/i] + (1 + cs — b) 



16 

: I 

7 1 



■ 47 



Qt 



b 

2~s 



) 2 

s xp 



It is clear that the coefficient of 3 pp [/i] can be made negative by choosing k 
sufficiently large, for all values of the other constants. To make the coefficients 
of 3 xx [h} 7 Q 2 , Q 2 negative we require that 



b > 1 + cs , b < 



f s + 4 7 



This is possible as soon as 
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a > (1 + cs)(ys + 47) 



and 



b < 2cs . 



c > s 



If we further require that a > As 2 c 7 then 2cs < \fac and therefore b < 2cs 
implies b < y/ac as well. This completes the proof of the proposition. □ 



5.3 Completion of the proof 

To complete the proof of Theorem [T] we appeal to Assumption [3] Using the 
logarithmic Sobolev inequality (|27|) in (|50|) we obtain 



^-£\h] < -(2da)S)[h] 
dt 

and combining with the second inequality in (|49[) we infer that there exists a 
constant A > such that 

p[h] < -XS[h] . 

Whence £[h] < £ [/i; n ] exp(— \t) < C(3 xx [h in ] + 3 pp [/i; n ]) exp(— Xt) and by the 
lower bound £[h) > k®[h], see Lemma [9j the entropy decays exponentially as 
well, which is the main claim of Theorem [TJ 
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A Appendix: Cauchy's problem 



In this appendix we discuss the global existence and uniqueness of solutions to 
the Cauchy problem for equation ([12")) in the case when the dimensions of the 
spaces AT and M. coincide, i.e., N = M . The following discussion is based on the 
methods introduced in |171 Ch. 5], except that we work in a different functions 
space. To adhere with the conventions used in [17], we rewrite (fT"2]) as 

d t h + Ah = , 

where 

A = -A p - W + v(p) ■ = —L + T . 

The domain D(A) of the operator A is chosen as the space of C°° functions 
on T N x WL N with compact support in the p € variable, which is dense in 
% := L 2 (dxdfi) . Our first purpose is to prove that the closure of A generates a 
contraction (dissipative) semigroup in %. To this end we need to assume that 
the quantities g, v, E are C°°. Furthermore we assume that 

-*■(), as \p\ oo = 1,...N. (59) 

We divide the proof in three steps. 

Step 1: A is accretive. 

By m, 

<h\Ah> n = - <h\Lh> u + <h\Th> u = J g(d p h 7 d p h) > . 
Recall that all integrals are extended over x ~R N with measure dxdjx. 



Step 2: A is hypoelliptic. 

Let a = \J g~ x (i.e., the positive definite matrix such that a 2 = g~ l ). A straight- 
forward calculation shows that the operator —A can be written in Hormander's 
form: 

N 
i=l 

where 

Yoh = (divp a) ■ aV p /i — g{d p E, d p h) — Th , 
Y^h — a^dpkh . 
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To prove that the operator A is hypoelliptic, we will show that —A satisfies a 
rank 2 Hormander's condition, namely that the vector fields 

Y (i) , z w :=[r ,y (i) ] 

form a basis of R 2N . To this purpose we observe that 

% = B*0p» + c/a^ = p (fc) + clx {1) , 



where 



C{ = a\d p ^) 



and B is a p— dependent N x N matrix, whose exact form is irrelevant for what 
follows. Thus the linear transformation {Xrj\, Pu)} — >• {Y(i), Z^} is represented 
by the matrix 

/ a 
\ C B 

The determinant of F is given by 

|detF| = deta|detC| = det g\ det{d pk v {1) ) | , 

which is positive because det g > and, by Assumption[2j the determinant of the 
matrix d p kv^ is non-zero. Thus {Yu\,Zu\} is a new basis of R 2N , concluding 
the proof. 

Step 3: The closure of A is maximally accretive. 

By [13 Th. 5.4] (see also [H]), it is enough to prove that the range of A + A is 
dense in H, for some A > 0. We need to show that if h € H is such that 

<h\(X + A)f> H =0, for all / e D(A) , (60) 

then h = 0. Note that implies that h is a distributional solution of 

(X-L-T)h = 0. 

Since the operator in the left hand side of the latter equation is hypoelliptic, 
then we may assume that h £ C°°. Let us begin by proving that the following 
identity holds: 

A J tfh 2 + J g{d p {(t>h),d p {ct>h)) = J e g {d p ct>,d p ct>)- J h 2 ^, (61) 

for all 4> E D(A). To prove (|MT) . we use that, by (i)-(ii) of Lemma [TJ 
(A + A)(hf 2 ) = /i(A + A)h + hAh - 2g(d P h,d p f 2 ) , for all f u f 2 e C°° . 
Setting fx = <p, f 2 = 4>h and multiplying by h the resulting identity we get 
4>h{\ + A){<j)h) = h(X + A){4> 2 h) - h 2 ^ + 2hg(d p (j), d p ((f>h)) . 
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Integrating and using that <h\(\ + A){4> 2 h)>u— 0, by ([SO]), we have 



/ <ph(X + A)(4>h) = J h 2 4>L4>~ J h 2 (j)T(l> + 2 J hg(d P 4>,d p (4>h)) . 

Using (fTU)l in the l.h.s. and in the first term in the r.h.s. of the previous identity 
completes the proof of (|61j) . Now let k g N and choose a family of test functions 
4>k of the form 

<Pk(x,p) = ^{p/k) , 

where i/j € C c °°, < ip < 1, V = 1 for p G 5(0,1/2) and suppV> C 5(0,1). 
Whence Ttfik = 0. Substituting in (TBTT) we obtain 

\J<t>lh 2 <^J h 2 g»d p ^d pi ^ X \p\<k ■ 

Having assumed (|55|) . we obtain 

A / 4> 2 h 2 < e(k) , 



where e(k) — > as fc — > oo. This finally entails that ft, = 0. 

We may now sketch the proof of the global well-posedeness of the Cauchy prob- 
lem in the setting of L 1 solutions. 

Theorem 3. Let g,v,E be C°° functions such that ([55)1 holds. Given < h ln g 
L 1 (T JV x R w ; dxdfi), there exists a unique h € C([0, oo); L 1 (T Ar x R^jdxd/i)) 
solution of (|10p toztft. initial datum h- m . 



Proof. Approximate the initial datum by a sequence /ii n)fn of smooth, non- 
negative functions in the domain D(A). By the preceding result, for each 
fixed m G N there exists a unique h m £ C([0, oo), L 2 (dx dfi)), solution of (TlU)) . 
Moreover by standard methods (see [TU1 [21] for instance) one can prove the 
^-contraction property: \\hk — ftmlU 1 < ||Mn,fe — ^hvrolU 1 - Thus the sequence 
/i TO converges in L 1 to a solution with the regularity stated in the theorem. 
The uniqueness is also a consequence of the ^-contraction property. The non- 
negativity of solutions can be proved by studying the evolution of a suitable 
regularization of sign(/i) (see again [TQl E3])- ^ 



B Appendix: Validity of the Logarithmic Sobolev 
inequality 

In this appendix we provide a sufficient condition for the validity of the loga- 
rithmic Sobolev inequality (|27jl . 
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Let Ajj denote the matrix inverse of A IJ = g(d p v( I \d p v( J ' > ), i.e., A IJ Ajk = 
5ft. We define the Riemannian metric G on M. x M as 

G = g lJ dp i <g) dp> + Audx 1 ® cfe J . (62) 

Moreover we define the vector field Q £ X(M) as 

Q = W-d p log Vdet(A/j). (63) 

Theorem 4. TTie inequality (j2"7| /loMs i/ 

Ric G (Z,Z) - (V G Q*)(Z,Z) > aG(Z,Z) , for all Z £ X(M x JV) , (64) 

where Ric G is i/ie Ricci tensor of G and V G is t/ie covariant differential asso- 
ciated with G. 

Proof. Consider the non-degenerate Fokker-Planck equation 

ft/ = A G / + Q/ (65) 

on M. xM . The entropy functional and entropy dissipation functional associated 
with (|65|) are exactly £> and 3 = Jpp+O^. Moreover ([64]) asserts that the metric 
G and the vector field Q verify the curvature Bakry-Emery bound condition. 
Thus the logarithmic Sobolev inequality (|2"T1) follows by the results in [6]. □ 



Using the relations between Ric G , V G and Ric,V p , the bound (|64"|) can be ex- 
pressed in terms of inequalities on the quantities g,W . These inequalities are 
in general very complicated, unless Au enjoys some simple structure, as in the 
statement of Corollary [1] below. 

Corollary 1. Let Ajj — the matrix inverse of A IJ = g(d p v^ I \d p v^) — be of 
the form 

Au = £(p) 2 6u , (66) 

for some smooth function £ : M — > (0, oo). // there exist two constants n\ > 
H2 > such that 



N 



Ric- y;\ogu- -^dpC^dpC 



A„logC 



C 2 

g(d p u,d p () 



(X,X)> Kl g(X,X), 



< k 2 , 



for all X,p£ R , i/iera ([64]) Wds. 



Proof. For the proof we observe that when Aj j has the form (|66l) , the Rieman- 
nian manifold (.M x AT, G) is the warped product of the manifolds (R , 5) and 
(T w ,(5), where <5 is the flat Euclidean metric on the torus. See [22j Ch. 7] for 
an introduction to the geometry of warped product manifolds. In particular, 
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by Corollary 43 of [22, Ch. 7] we have that, for all horizontal (i.e., tangent to 
M) vector fields X, Y and vertical (i.e. tangent to J\f) vector fields V, W, the 
following identities hold: 

Ric G (X, Y) = Ric(X, Y) - jV 2 p ((X, Y) , 
Ric G {X, V) = 0, 

Ric G (V, W) = - + {N- l)^|lM)J G(V, W) . 



Moreover the vector field (|63[) becomes 

Q = d p {\ogu-N\ogO 

and by Proposition 35 of [351 Ch. 7] we have 

V G Q*(X,Y) = V P Q4X,Y) = (V 2 p logu-NV 2 p log() (X,Y) , 
V G Q*(X, V)=0, 

V o Q4V} w) = m G(Vj w) = ( 9(d P u d p Q _ N 9(d P C ; d P Q \ Qfy w) 

C \ (u Q J 

Now let X = Y = Zi M , the projection of Z onto Ai, and V = W — Z\ M , the 
projection of Z onto TV. Writing Z = Z\ M + Z\ M we have 

(Ric G - \7 G Q*)(Z,Z) = (mc- jV 2 p ^J (Z lM ,Z lM ) 

^ + (n-i) 9 -^^)g(z w ,z w) 



(V>gu-tfV>gC) (Z lM ,Z lM ) 

g{d p u, d p Q AT g(d p C,d P Q \ ^ /r7 . 
N ^ ) G i Z U> Z U) 



N 



Ric-V^log U --a p C®^Cj (Z\ M ,Z ]M ) 
A P logC + ^^) G (Z k ,Z k ). 

Under the given assumptions we have 

(Ric G — V G Q»)(Z, Z) > Ki G(Z,Z)-k 2 G(Z w ,Z w ) 

= ( Kl - k 2 )G(Z, Z) + k 2 G(Z 1m ,Z lM )>( Kl - k 2 )G(Z, Z) 

and the conclusion of Corollary Q] follows. 

□ 
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An example that is covered by Corollary Q] is the Fokker-Planck equation ^ 
with the classical velocity field v(p) — p and an isotropic diffusion matrix, i.e. 

D ij (p) = Il( P ) 2 5 ij , 

where II is a positive function. In particular g 1 ^ = II(f>) 2 £ y . Thus, since = 
p 1 (all indexes run from 1 to N in this example), we have 

A IJ = gftdptvWdpjvW = g ij Sl5l = g IJ = n(p) 2 £ /J , 

and so §^ holds with C(p) = l/n(p). 
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